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ABSTRACT 

Using the gapless second-order theory, the excitation energies and 
dynamics of the collective excitation for a partially condensed 
harmonic trapped quasi two dimensional bosonic gas has been 
studied. Our evaluated result for temperature dependent mean- 
energies show that the Kohn theorem is quite accurately satisfied 
for temperature T<0.8 T c . 
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INTRODUCTION 

The partially Bose-Einstein 
condensate trapped atomic gases provide an 
excellent test bench for developing finite 
temperature quantum theories. These weakly 
interacting systems can be modeled from 
first principles, and the experiments yield 
accurate and detailed information for 
comparison. Especially, the energies and 
decay rates of low-energy collective 
excitations have been measured tests for 
theoretical models. 

For dilute condensates at 
temperature much lower than the condensate 
temperature T c ' the Bogoliubov approxi- 
mation consisting of the Gross-Pitaevskii 



(GP) equation for the condensate wave 
function and the Bogoliubov equations for 
the quasi-particle excitations has proven to 
be accurate in describing the collective 
modes of the system. For higher 
temperatures one has to take into account the 
effects of the thermal gas component. 
Developing a theory that is computationally 
feasible and correctly models the system at 
temperature approaching T c is a challenging 
task. The most commonly used finite- 
temperature theory is the Hartree-Fock- 
Bogoliubov-Popov (HFB-Popov) approxi- 
mation. It neglects the dynamics of the 
thermal gas and the modifications in particle 
correlation induced by the condensate, but 
predicts quasi-particle energies in fair 
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agreement with the experiment. 1 The energy 
of the quadruple modes having azimuthal 
angular momentum quantum number qe=±2 
deviates from the theoretical prediction for 
temperature above 0.6T C , but lately this 
deviation has been interpreted to mainly 
arise from improper modeling of the time 
dependent external potential used in the 
experiments to excite the collective modes. 2 
In order to take into account the leading 
order quasi-particle interactions and the 
correlation induced by the condensate in the 
inhomogeneous case, several theoretical 
approaches have been suggested. 310 The 
dynamics of the condensate and the thermal 
gas has also been studied using various 
kinetic theories. 11 ' 15 The second-order theory 
for inhomogeneous, partially condensate 
gases presented in refs. 9 ' 10 uses systematic 
perturbation theory to take into account the 
interaction terms in the Hamiltonian. 
Recently, this theory was extended to take 
into account the time -dependent external 
perturbation used to drive the system in the 
experiments, leading to an agreement with 
the measured energies and the damping rates 



of the collective modes. 2 ' 16 

The second-order theories are 
computationally challenging, and there have 
been only a few numerical investigations of 
their prediction. 2 ' 17 ' 18 In this paper, we 
calculate the spectral distributions of the 
quasi-particle energies for a partially 
condensed Bose-Einstein Condensate (BEC) 
and compare the quasi-particle energies to 
the HFB -Popov results as function of 
temperature. Especially, we analyze the 
quasi-particle dynamics implied by the 
spectral distributions, observing that some 
collective modes should exhibit notable 
collapse and revival effects in trapped 



condensates. The possible existence of this 
phenomenon has been pointed out 
previously in Ref. 9 ' 10 , but it has not been 
studied in detail before. The collapse and 
revival of the excitations indicate that the 
energies and the damping rates alone do not 
suffice to describe the dynamics of these 
modes, i.e., the commonly used damped 
sinusoidal fit to the experimental data may 
be sufficient to describe the longer term 
dynamics of some modes. 19 24 

2.0 Mathematical formulae used in the 
evaluation 

Second-Order Theory 

In this section, one presents the 
second-order formalism for calculating the 
quasi-particle spectral distributions for a 
partially condensed, dilute, trapped BEC at 
finite temperatures. The starting point is the 
usual second-quantized Hamiltonian for 
structure less bosons 



** = H{H SP \\ J ) a l a ) \y\km)a}a)a k a m 



Where the creation and annihilation 
operators for a particle in state li) are denotes 
by a and respectively. The single particle 
Hamiltonian is given by the sum of the 
Kinetic energy and the external trapping 
potential as 



H SP = 



h 2 v 2 

2m 



+ V (r) 



and the dominant s-wave scattering at low 
temperature can be modeled by the effective 
low energy interaction potential 
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V(r) 



47Tah 2 S(r) 



m 



ijkm^O ^ 



V 



km)a]a]a k a m 



where a is the scattering length and m the 
atomic mass. This effective potential is 
inapplicable at high energies and leads to 
ultraviolet divergences in the theory which 
have to be renormalized in a proper way. 

We choose to use canonical 
ensemble with fixed total number of 
particles N. By defining the bosonic number 
conserving operators 



and S = N-(N 



is the number 

fluctuation operator of the non-condensate 
particles. The symmetries elements of the 
two-particle interaction potential V(r) are 
defined as 



ij |v * | km) = — (ij \V | krnj + (ij |v| krnj 



a 



(#o+l)' 



-1/2 



a,, 



where the index 0 refers to the condensate 
state and N 0 = dla 0 , one can write the 
Hamiltonian (1) as 



^o=I^+0 
Where, 



/=o 















V 







H 0 = N Q 



H 



SP 



H, =N, 



|0> + |iV 0 (0|V^|0 



\\0) + N 0 (i0\v s \00) a] + H.c. 



**2 = Z [i (H SP I j) ~ Hj + 2N 0 (Oi \V S I JO) a)a i ] 



ij*0 



Hij\V'\00)a]a)+H£. 



# 3 = Z 4N^(ij\v s \K0)a*a]a k +H.c. 

ijk*0 



and as 



0 



H 



j) + N 0 00 V 100 



Where the average number of atoms in the 
condensate state is given by 



N 0 = N-{N ex j . Above the average {...} 

refer to quantum expectation values and H.c. 
for Hermitian conjugate. 

In the zeroth-order approximation, 

one solves the ground state |0) of H 0 

alone, which makes the linear Hamiltonian 

H { \ vanish. The excitations are found in 

lowest order by diagonalizing H 2 and the 
number of the condensed particle N 0 has to 
be tuned such that the total number of 
particles satisfies N = N 0 + N ex , 

It is convenient the use an 
orthonormal single -particle basis 

£ii r ) = \ r \i) f° r a ^ * ~ 0' wnere 
£ 0 (r) is the condensate wave function 
given by the Gross-Pitaevskii equation 
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VX 0 (r)=V„Jr)£ 0 + N 0 U 0 \£ 0 (r)\ 2 £ 2 =A£ 0 (r) p ( r ) = ^U u . (rtf +\v, (r)| 2 ~|n ( . (r)| 2 | 



with f/ 0 = A7Tafi 2 1 m. . The GP equation is 

obtained by rninirnizing w i tn res P ect 

to £ 0 (r). Diagonalizing H 2 using the 
Bogoliubov transformation 



' L(r) 


M(r)V M; (r) 




V -M *(r) 


-L(r)J[v,.(r) 





where 
and 

«i('-) = Z Al oW(r) 

are the quasi particles. G ; the quasi-particle 
energies, and operator 



*(r) = 2« < (r)»;(r)(2n l +l) 

One finds the perturbative Hamiltonian 
H pert =AH 0 +AH 1 + AH 2 + H 3 + H 4 

Where the nonquadratic term H 3 and H 4 
are to be calculated using the improved 
condensate wave function. Note that our 
notation AH ; differs somewhat from that in 
Refs. 9 and 10. 

The perturbation term AH 0 is just a 
real number and can be easily taken into 
account. In addition to it, in first-order 
perturbation theory only the terms AH 2 and 
H4 containing even numbers of quasi- 
particle operators contribute to the energy 
shift. 



E pert (s,l) = 



H 



per! 



L(r) = H^-A + 2N 0 U Q \C-(r)\ 2 



where \s) is a quasi-particle occupation 

number eigenstate. In second-order 
perturbation theory, one can in fact neglect 
the terms AH 2 and H 4 because it turns out 
and M (r) = NJJ a Q (r) have been that their contribution is of the same order as 

introduced. The quasi-particle amplitudes 
must satisfy the orthogonality and symmetry 
relations 



the contribution of the other terms in third- 



order perturbation theory, 
needs to calculate only 



Thus, one 



Jdr[i* < (r)i*;(r)-» < (r)tt*;(r)] = < y tf 
\dr\u t (r)v j (r)-V i (r)vUj(r)~\ = 0 



r Atfj + tfj 



E-E. 



The quasi-particle energies are 
We express the particle operators in terms of calculated as total energy changes in the 



the quasi-particle operators, yielding 



system which the corresponding quasi- 
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particle occupation number by 1, while the 
total number of particles is held constant. 
This yields the corrected excitation energy 

E p (Z') =g p +AE: + AE> + AE{ + AE 3 P (?) 



. , -z 2 /2(2fl. 2 ) / 77J 

where <7(z)e I \ja z 7T is a 

Gaussian profile and a t = I mCD i are the 
harmonic oscillator lengths of the trap. 

Table Tj 

An Evaluated results of temperature dependent 
mean energies of the excitation modes using 
second order theory and also with HFB-Popov 
(Hartree-Fock-Bogoliubov) theory. The exact 
energy hco r of the Kohn modes are also given 



where the A terms are given in Eqs. and the 
complex energy parameter z' should not be 
mixed with the fugacity. Calculating the 
excitation energies as function of z' yields 
the dynamics of the excitations in the 
following way. The time evolution operator 
U(t) of the system may be written in terms 
of the Fourier transform of the resolvent 
operator 

g{z') = (z'-h)~ 1 
h 



U(t) = --t oo e iC0 'Im 6{hco-i0) 
NUMERICAL METHODS 



da 



One considers a pancake-shaped 
system in a harmonic potential 

Vtrap {r) = \rna>y + ^mO) 2 y 2 + ^m(0 2 z 2 

where the trapping frequencies are 
G) r = Q) x = Q) y and 0) z with G) z » 0) r . For 

sufficiently strong trapping potential in the z 
direction, the condensate wave function and 
the thermo-dynamically relevant quasi- 
particle amplitude can be approximated to be 
in cylindrical coordinates (r, 8, z) of the 
factorized form 

and 

u{r) = u{r)a{z)e (q ° +m)e 
V{r) = v{r)a{z)e l( - qe+mf 





E/hco r 


T/Tc 


Second- 


HFB-Popov 


Energy of 




order 


Theory 


Kohn 




Theory 




Theory 


0.0 


2.058 


2.153 


1.079 


0.1 


2.092 


2.224 


1.082 


0.2 


2.145 


2.278 


1.097 


0.3 


2.181 


2.324 


1.116 


0.4 


2.226 


2.359 


1.132 


0.5 


2.274 


2.396 


1.138 


0.6 


2.345 


2.438 


1.142 


0.7 


2.378 


2.473 


1.145 


0.8 


2.455 


2.524 


1.154 


0.9 


2.481 


2.545 


1.166 


1.0 


2.550 


2.607 


1.178 



DISCUSSION OF RESULTS 

In this paper, we have evaluated the 
temperature dependent mean energies of 
excitation of Bose-Einstein Condensate. For 
numerical estimation, we have taken a model 
of pancake - shaped cloud consisting of N = 
2000 23 Na atoms trapped with trapping 
frequency ro c = 25 x 350 Hz. The radial 
trapping frequency co r - m x = ro y may be 
chosen freely with only the constraint ro z > 
co r . The parameters are chosen from ref. 23. 
We have taken the theoretical formalism of 
M. Mottene, S. M. M. Virtaken and M. M. 



Journal of Pure Applied and Industrial Physics Vol.2, Issue 2, 1 April, 2012, Pages (84-141) 



131 Md. Shabbir, et al, J. Pure Appl. & Ind. Phys. Vol.2 (2), 126-131 (2012) 



Salomara. In table Ti , we have shown the 
evaluated results of temperature dependent 
mean energies of the excitation using second 
order theory together with HFB -Popov 
theory and exact energy of Kohn modes. 
Kohn modes are also called center of mass 
oscillation modes. 26 According to Kohn 
thorem 27 a system of harmonically trapped 
interacting particle in any eigenstate of the 
Hamiltonoian has an eigenstate with the 
amount hcoi . The Bogoliubov theory, in 
which the thermal gas components is 
neglected, implies Kohn modes to have this 
exact energy. In higher order theories, the 
dynamics of the thermal gas and its 
interaction with the condensate have to be 
taken into account accurately obtained 
results in agreement with the Kohn theorem. 
In our calculation, it is shown that within the 
second order theory, the energy of the Kohn 
mode is very close to hro r for temperature T 
< 0.8 Tc. 
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